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We study the critical behavior of a random field 0{N) spin model with a second-rank random 
anisotropy term in spatial dimensions 4 < ci < 6, by means of the replica method and the 
expansion. We obtain a replica-symmetric solution of the saddle-point equation, and we find the 
phase transition obeying dimensional reduction. We study the stability of the replica-symmetric 
saddle point against the fluctuation induced by the second-rank random anisotropy. We show that 
the eigenvalue of the Hessian at the replica-symmetric saddle point is strictly positive. Therefore, 
this saddle point is stable and the dimensional reduction holds in the expansion. To check 
the consistency with the functional renormalization group method, we obtain all fixed points of the 
renormalization group in the large A'' limit and discuss their stability. We find that the analytic 
fixed point yielding the dimensional reduction is practically singly unstable in a coupling constant 
space of the given model with large A'^. Thus, we conclude that the dimensional reduction holds for 
sufficiently large A'^. 
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I. INTRODUCTION 

The random field O(A^) spin model is one of the sim- 
plest models with both a site randomness and a short 
range spin correlationi Despite intensive research for 
about three decades, our understanding of this model 
is not yet satisfactory (for recent review, see Ref. f3)- 

Dimensional reduction^ is one key to clarify the na- 
ture of this model. Dimensional reduction claims that 
the critical behavior of the d-dimensional random field 
O(A^) spin model is the same as of the (d— 2)-dimensional 
pure 0{N) spin model, where d is the spatial dimension. 
It has been shown by rigorous proofs^^ and numerical 
calculations of critical exponentsSi^iS that the predic- 
tion of dimensional reduction is incorrect in the random 
field Ising model below four dimensions. In dimensions 
more than 4, however, the critical phenomena of the ran- 
dom field 0{N) spin model should be further studied. 
In particular, the breakdown of the dimensional reduc- 
tion and the possibility of an intermediate phase between 
the paramagnet and ferromagnet phasess are still under 
controversy. 

Mezard and Young considered the possibility of the 
glassy phase by replica symmetry breaking. They dealt 
with the random field 0{N) 0'* model, and studied the 
critical behavior by using the replica method and the self- 
consistent screening approximation (SCSA), which is a 
truncated Schwinger-Dyson equation for a two-point cor- 
relation function. Under the assumption of replica sym- 



metry, the dimensional reduction appears and the critical 
exponents of the connected and disconnected correlation 
functions rj and fj satisfy fj = rj. They showed that the 
replica-symmetric correlation function was, however, un- 
stable as a solution of the SCSA equation at 0(l/iV). 
They proposed a replica-symmetry-breaking correlation 
function, where they found 2ry ^ fj. Following Mezard 
and Young, the instability of the replica-symmetric solu- 
tion against replica symmetry breaking has been reported 
in several papers^i*i^ However, the physical meaning of 
the instability in the SCSA equation is still unclear. 

Fisher and Feldman pointed out the breakdown of 
the dimensional reduction due to the appearance of 
the infinite number of relevant operators near four 
dimensions. ^■^'^^ Fisher showed that all possible higher- 
rank random anisotropics are generated by the func- 
tional renormalization group recursion relations for the 
0{N) nonlinear a model including only the random field 
term. The random field and the random anisotropics 
are marginal operators in d = 4. Then he treated the 
nonlinear a model with a random field and all the ran- 
dom anisotropy terms, and calculated the one-loop beta 
function for a linear combination of them in d = 4 -I- e 
under the assumption of replica symmetry. He showed 
that there is no singly unstable fixed point of 0(e) which 
gives the results of dimensional reduction, and that the 
flow goes into the regime where nonperturbative effects 
are important. Therefore, he concluded that the dimen- 
sional reduction breaks down at least near four dimen- 
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sions. Feldman carefully reexamined the one-loop beta 
function obtained by Fisher. He treated a differential 
equation as the fixed point condition and found nonan- 
alytic fixed points which control the critical phenomena 
instead of the analytic fixed ones. He calculated the ex- 
ponents 77 and 77 for = 3, 4, 5 in 4 + e dimensions nu- 
merically; then he concluded that dimensional reduction 
breaks down near four dimensions for several finite N. 

These studies indicate the breakdown of dimensional 
reduction in the random field 0{N) spin system. How- 
ever, the relation between the renormalization group and 
simple 1 /A^-expansion methods has never been discussed. 
Thus, it is important to study the relation between the 
stability of the replica-symmetric saddle point and the 
analytic fixed point in the functional renormalization 
group for large N. 

In this paper, we study the random field 0{N) 
spin model including random anisotropy by a simple 
I /N expansion and the functional renormalization group 
method. We study the robustness or fragility of the sys- 
tem against the random anisotropy perturbation. First, 
we study the stability of the replica-symmetric saddle 
point in spatial dimensions 4 < d < 6 by the simple 1/N 
expansion. To investigate the stability of the replica- 
symmetric saddle point against a small perturbation of 
the second-rank random anisotropy, we employ the crite- 
rion for stability used by de Almeida and Thouless.«iSi We 
find that the eigenvalues of the Hessian are strictly posi- 
tive, and the replica-symmetric saddle point remains sta- 
ble against the second-rank random anisotropy. There- 
fore, the dimensional reduction works well for large A^. 
Next, we check the consistency of this result with the 
functional renormalization group analysis in 4 + e dimen- 
sions. We solve the fixed point condition of the renormal- 
ization group in the large N limit, and study the stability 
of all fixed points. We solve the eigenvalue equation for 
the infinitesimal deviation from the fixed points. We find 
that the analytic fixed point yielding dimensional reduc- 
tion is singly unstable. Careful analysis of the eigenvalue 
equation for the infinitesimal deviation from this fixed 
point is done in terms of 1/A^ expansion. We find in- 
finitely many unphysical modes which should be elimi- 
nated. In practice, the analytic fixed point yielding the 
dimensional reduction is singly unstable for sufficiently 
large A'^. Therefore, our simple 1/N expansion is consis- 
tent with the functional renormalization group method 
and we conclude the dimensional reduction 

^^^^ N -2 

for sufficiently large A^. 

This paper is organized as follows. In Sec. ^1 we briefly 
review the large A^ behavior of the random field 0{N) 
spin model in the absence of random anisotropy. In Sec. 
mil we introduce the second-rank random anisotropy 
term, and perform the 1/N expansion for the random 
field O(A^) spin model with the second-rank random 
anisotropy term. We should integrate over the "off- 



diagonal" fluctuation introduced through the Hubbard- 
Stratonovich transformation for the second-rank random 
anisotropy. Solving the saddle-point equations under the 
assumption of replica symmetry, we have two solutions. 
Then we calculate the free energy densities at high tem- 
peratures in both solutions, and compare with the result 
of the high temperature expansion without the replica 
method. As a result, the solution is uniquely determined. 
Details of the calculations of the free energy at high tem- 
peratures without the replica method are relegated to 
Appendix We also calculate the critical line, and the 
eigenvalue of the Hessian at high temperature and near 
the critical point. The stability of the replica-symmetric 
saddle point is investigated. Details of the calculations of 
the eigenvalue of the Hessian are relegated to Appendix 
IbI In Sec. IIVI we compare our results with those of 
a renormalization group study. Technical details of the 
renormalization group for large A'^ are presented in Ap- 
pendix El Finally in Sec. we summarize the results 
obtained in this paper, give some comments on the crit- 
ical phenomena of both the lower and the upper critical 
dimensions on the basis of the results, and mention fu- 
ture problems. Calculation of loop integrals is exhibited 
in Appendix IdI 

II. CRITICAL BEHAVIOR OF RANDOM FIELD 
0{N) SPIN MODEL IN THE LARGE TV LIMIT 

In this section, we briefly review the 1/N expansion for 
the 0(A'^) spin model with only a random magnetic field 
under the assumption of replica symmetry. The stability 
of the replica symmetric saddle point is studied. 

We consider the random field O(A^) spin model on a 
d-dimensional hypercubic lattice with the lattice spacing 
unity. Let L be the linear length of the d-dimensional 
hypercubic lattice, and V the number of lattice sites {V = 
L'*). The Hamiltonian is given by 

H = -jJ2Scc-Sy-J2hx-S^. (I) 

Here (x, y) denotes the summation over the nearest 
neighbor pairs of the lattice sites x and y. J is the ex- 
change interaction, and we take J > 0. Sx denotes an 
A'^-component spin variable on the site x with a fixed- 
length constraint = 1, and denotes a Gaussian 
random field with zero average. Taking the average over 
the random fields {h^} by using the replica method, we 
have the following replica partition function: 

X /"(n n ^dSx,^6iSx,J - I)]e-^^-, (2) 

^ X Q = l ^ 



X a, 13 
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Here A^; stands for the lattice Laplacian. In the momen- 
tum representation, the lattice Laplacian is represented 
by Afe — 2^^j^^(cosfc^ ~ !)• /3 is the inverse tempera- 
ture, and T is the temperature; /? = 1/kT. A denotes 
the strength of the Gaussian random field. The replica 
index is denoted by a = 1, n. We rewrite S{Sx,a^ — 1) 
in terms of the auxiliary variable Xax- 



SiSx,a'^ - 1) = 



An 



-l)/2 



(4) 



After integrating over the spin variables {Sx,a}i the 
replica partition function becomes 



nV /r,_Ar\ NnV/2 

\4Tr J \ P ) 



-Sctf 



(5) 



Ses = -^(a:|Trln(-JA,l„ + ;c)b 



where 1„ is an n x n unit matrix, and x is n x n 
symmetric matrix with 



Xaflx — iXaxSaf3 — /^A. 



(7) 



We study the large N limit below. The large N limit 
is taken with NT (or jS/N) and NA finite. Then, we 
redefine the parameters as follows: 



7VT->T (^^^P 
NA A. 



(8) 



A. Saddle-point equation and replica-symmetric 
approximation 

Differentiating SoS by iXax, we get the saddle-point 
equation 



SS,s N 



SiXo 



1 



Here we assume the replica symmetry 
iXax = . 

In this assumption, 

) 



.)-f ^O.(ll) 



(12) 



1 



-JAk + m- 



f3A 



-JAk + m? 



^G§,<5a/3 + /3AG[J,^G^f. 
The saddle-point equation becomes 

1 = ^a{m^) + Ab{m^), 



(6) where 



a(m^) = — — . 

V ^ - JAk + m^ 



i-JAk + m^)^ 



(13) 
(14) 

(15) 
(16) 



In the thermodynamic limit V — > oo, a{m?) and h{m?) 
change over the integrals: 



aim ) = 



2 \ ^ — 



d'^k 1 

d'^k 1 



kei-.,.]" (27r)'^ (-JAfe + m2)2 



(17) 
(18) 



Thus, the replica partition function is rewritten as fol- 
lows: 



Z = e 



{4n J {f3 



TV 



Scs = y ^(a;|Trln(-JA,l„ + x)b 

X 



(9) 



(10) 



d'^k TT f r dk^ 



fce[-7r,Tr]'' (Ztt)*^ \J-TT 



2tt 



(19) 



Near the critical point, m becomes small, and then the 
integrals H17() and 118() can be expanded in terms of m 
for 4 < d < 6, 



(20) 
(21) 



where ao, ai, bo, and 6i are positive constants. The 
derivation of these is shown in Appendix D. Inserting the 
above expansions into the right hand side of Eq. H14|l . we 
get 



1 = kT{ao — flim^) + A(feo ~ bim' 



d-4\ 



(22) 
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At first, we study two special cases: one is the A = 
case, and the other is the T = case. Putting A = 0, we 
have the following expression for the saddle point to: 



TO 



flo T - Tf 



(pure) 



ai 

^y(purc) - 



T 



1 



This indicates 

V 



TO 
1 

2' 



(T - rjp 



urc) \—v 



(23) 
(24) 

(25) 
(26) 



This result is identical with that of the mean field theory 
of the pure system as expected. In the case of T = 0, to 
is expressed as follows: 



feo A - A 



(T=0) 



6i 

. (T=0) ^ J_ 
6o 



This indicates 

J/ = 



(A_Af=o))-, 



1 



d- 4 



(27) 
(28) 

(29) 
(30) 



In d = 4 + e, the above result agrees with that of pure 
systems in d = 2 + e in the leading order. Next, we study 
the case of T 7^ and A 7^ 0. The saddle point to is 
expressed as follows: 



TO 



1 

A6l 



(fcTao + A60 - 1; 



l/(rf-4) 



(31) 



Putting TO = 0, we can get the critical line between fer- 
romagnetic and paramagnetic phases: 



, (T=0) 



(pure) 



= 1. 



(32) 
(33) 



The phase diagram is depicted in Fig. ^ to is rewritten 
by using Tc and Ac as follows: 



TO 



This indicates 

V 



1 1/(^-4) 



A61 ^j-^Cpui-c) aF"°^ 
(T-T,)i/('i-4) (A = A,), 
(A- Ac)i/('^-4) (T^Tc). 



TO-^ ^ (A- A,)- 



(T-Tc 



d-4 



(34) 

(35) 
(36) 



In d = 4 + e, this result is identical with that of pure 
systems in d = 2 + e in the leading order. 




FIG. 1: Phase diagram of the random field 0{N) model 

B. Stability of replica-symmetric saddle point 

We put 

XafSx = (to^<5q/3 - PA) + ieaxSa/S 



Xafi + ^Xafix-, 



(37) 



and expand the effective action ^cff up to the second 
order of ^Xafix- The second-order term of ixa^x for the 
effective action S'off becomes 



Tr- 



:Sx 



1 



-JA^ln+X -JA^ln+X 



Sx 



N 



4 Jke\-Tv,TT]'i (27r) 



N 



a,f3 



-ki-i-kt-k, 



4 Aehrr,^]^ (27r)'^ 
in the thermodynamic limit. Ilapk is 



(38) 



n 



aPk 



(27r)'' 0'=-' °' 
= [{A * A)k + {A * B)k + [B * A)k]5o.p + {B * B)k 
= ncfc(5a/3 + ndfc, (39) 



(A * A)k 
{A*B)k 



d'^a 
{B * A)k 

d'^q 



(3A 



[-Tr,TT]<i (271') 



(40) 



(41) 



{B*B)k - (/3A)2 / ^Gg, (42) 
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The expression is an n-dimensional vector whose ele- 
ments are eak- 



(43) 



and life denotes an n x n matrix whose elements are IIq 



(44) 



Calculating the eigenvalues of the matrix 11^, we have 



(45) 



Taking the n — s- limit, we can obtain the following 
expression for the eigenvalue: 



lim Afe „ = ricfc 



h7r,7r]'i (271") 



2/?AGoj._ Gg ) 



(46) 



Therefore the eigenvalue Xk is positive for T > and all 
A:. This result indicates that the replica-symmetric saddle 
point is stable against "diagonal" fluctuations Cak, and 
therefore it is possible to integrate out the fluctuations 

As seen in Eq. Ij^Tfl . the field includes no "off- 
diagonal" terms. In the next section, we shall study the 
effects of a second-rank random anisotropy on the crit- 
ical phenomena of the random field 0{N) spin model. 
We will find that the off-diagonal fluctuation is intro- 
duced through the Hubbard- Stratonovich transformation 
for the second-rank random anisotropy term. 



C. Calculation of r/ and r; in 4 + e dimensions 

Here, we calculate the critical exponents i] and fj. For 
simplicity, we put J = 1. At criticality, the lattice Lapla- 
cian becomes Aj, — —k^. Eq. (|39|l is 

Ha^fc = (co + cik"-^ + C2k''-'')5^fi + c^k''-\ (47) 



for m? = 0, where 



Co 



ci = 



1 r(^)[2r(i^^)2-ir(^)2] 



(47r)'i/2 



r(d-4) 



C2 



C3 



2A T{^) ^(d-2\^(d~A 



(47r)''/2 r(d - 3) V 2 



r 



A2 r(^)^^d-4^' 



(4^)d/2 r(d-4) 



(48) 



Let us compute the correlation function at the second 
order of the perturbation. Up to the second order of 
tak, we get the following expression for the correlation 
function: 



^ a — 1 



X exp(-(5^S'off) 



1 



a/3 



Op j (2jj-)d^0p-q\^iq^^,-qh'-^t 

7,(5 



■50 
Op 



' ^Op "^op '-'Op 



7,(5 



— ^Op 2-^ ^Op ^p ^p ' 
7,(5 

where Z^, and {e^qes-q)t are defined by 

Ze^(ji{ 2?ea) exp(-<525eff), 

y7(5 ^ /" J^niS I \ 
p ~ J ^27r)'' P^'' 



(49) 



(50) 



{^'yq^S, — q)e - 



z, 

2 

]VV Co +ci 9-^-4 + C2g'^-6 



e^ge5,-gexp(-J^S'off) 



(51) 



(co + cig''-4 + c2(7''-6)2^ 
In 4 < d < 6 and in low momentum, {£-/qes.-q)e becomes 



2 / 1 



\^iq^S-q)e 



6-d A 



Nc2\q'^~^ 2 g'' 

Thus, we get the following vertex function: 



— • (52) 



ip' + mi)S^^ - Ar + Sf - S;"^ (53) 
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where Sq"'', jtt-r, and Ar are defined by 



A'(6-d) 



m (5q/3 - A = m^daf3 - Ar, - Ej 



'a/3 



(54) 



Sq"'' does not include an infrared divergence. At the 
criticahty to^ = 0, we have 

(G;i)„0=p2,5„/5-AR + Ef -S^"^. (55) 
Ep'^ — Eq"'' is calculated as follows: 



NC2 

where Dp and i?p are 

A 1 d-6 
P - 2" (47r)'^/2 r(d/2) 



(56) 



2 1 2 

p Inp , 



A2 1 d-6 

" " ~ (4^)'*/2 r(d/2) 



Inp^ 



(57) 



Thus, we get the following expression for the vertex func- 
tion: 

{G-')a0 - [p' + J^^Dp^^p - (^Ar + -^^Epy{58) 
At criticality T — Tc, the correlation function behaves as 

1 „ Ar 



G 



af) 



(59) 



at low momentum; namely, the vertex function behaves 



r(-2?7 



~/(l - ?7lnp)(5a;3 - Ar[1 + (f; - 277) Inp]. 

(60) 

From Eqs. H58|l and (|60|l . we see that ry and 77 are of the 
order of in d = 4 + e as follows: 



77 = r; 



N 



(61) 



This result of rj is consistent with that of a pure system 
in d = 2 + e up to order e. The result fj — rj confirms the 
dimensional reduction. 



III. CRITICAL BEHAVIOR OF RANDOM 
FIELD 0(7V) SPIN MODEL WITH 
SECOND-RANK RANDOM ANISOTROPY IN 
THE LARGE- TV LIMIT 



anisotropy: 



rep 



^^5x.a(-JA,5„^-/3A)5. 



x,f3 



(62) 



The second term of the right hand side in the Hamil- 
tonian is the second-rank random anisotropy term, and 
g denotes the strength of the random anisotropy. The 
second-rank random anisotropy term is decomposed into 
diagonal and off-diagonal parts: 



~2~ '^^i^x,a ■ Sx.JlY 

a.J3 



(63) 



We rewrite the g term in terms of the auxiliary variable 
Qa^x as follows: 



X exp 



(64) 



We should note that the off-diagonal variable Qafix is 
introduced through the above transformation. Using the 
above equation and Eq. the Hamiltonian becomes 



/3^^rcp — 77 Sax{-J Ax^ap + Xa0x)S^ 



i3x 



X a,/3 



n n 

EE^^^ + i^EE^'/^- (65) 



X a—1 



X a=tfj 



where Xaf}x is 



In this section, we study the large N behavior of the 
following Hamiltonian including the second-rank random 



XaPx 



i^ax — PA (a = /?), 
-/3A-Q„;3. («^/3). 



(66) 
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After integrating over the spin variables {Sax}, the 
rephca partition function becomes 



Z 



An) [ P ) [ATrgJ 



^^nV(3^g/2^nV(3Jd 



(67) 



'S'eff = 



N 



2 J2{x\TTlni-JAxln + x)\x) 

X 

-fEE*^"^ + i:EE^?^/^^- (68) 

X a=l ^ x a^p 

The expression 1„ is the n x n unit matrix, and x is 
the n X n symmetric matrix whose elements are (|66ll . 
As in the previous section, we study the large N limit. 
The large N limit is taken with NT (or f3/N), NA and 
Ng staying finite. Then we redefine the parameters as 
follows: 



NA A, 
Ng^g. 



(69) 



Thus, the replica partition function is rewritten as fol- 
lows: 



Z = 



[attJ [(3) [iTTg) 

^^NnV0^g/2^NnV[3Jd 



n(Ti-l)V/4 



-Sotf 



(70) 



N s-^ 

Scs - -^(a;|Trln(~JA,l„ + x)|x) 



^EE^A"^ + SEE^?^..- (71) 



X a^l3 



A. Saddle-point equations and replica-symmetric 
approximation 

Differentiating ScS by iXax and Qaf3x respectively, we 
have the saddle-point equations 



cff 



N 



SiXax 



1 



Nj3 



(72) 



5Qa 



N 



fix 



- JAj.l„ + X/ aP 



N 



0. 



Here we assume the replica symmetry 

iXa — rm? , 
Qap — Q, 

(to2 + Q)S^^ -(3A-Q 



In this assumption, 
k 



M^S^p - ipA + Q)= xc.13. 



Qaf3x 

(73) 

(74) 
(75) 

(76) 



I 



-JAfel„ + X/ a/3 
1 



/3A + Q 



-JAk + AP^"^^ {-JAk + A'Py 
^G'^t,5ap + {pA + Q)G^^. 

The saddle-point equations become 



1 



a(M^) 



A + iQ)6(M2^ 



2g{l3A + Q)b{AP), 



(77) 

(78) 
(79) 



where 



a(M2 



V^oo 



-E' 



-JAk + AP 
d'^k I 



(80) 



(81) 



/fcg[_^,^]d (27r)'^ (_ JA;, + M2)2 
Thus, the saddle-point equations are rewritten as follows: 

(82) 
(83) 



= 2(ig[\~-a{M^)], 



Q = 2/35 



AbjAP) 
1 - 2gb{AP) ' 



We look for the intersections of these saddle-point equa- 
tions. For convenience, we define p = Q/{2f3g). Then, 
the saddle-point equations are rewritten as follows: 



p = I- ^a{M% 
Ab{AP) 



1 - 2gb{AP) 



(84) 
(85) 



The graphs of Eqs. (p^ and are drawn in Fig. El We 
find that there are two intersections at high temperature: 



Ab{Mt*) 
l~2gb(Mf) 

q = pX 



{AP = M2*), 
(M^^A/f), 



(86) 
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Ta{lvf) 



M 



/r2 



1 - 2gb{M^) 



FIG. 2: The graphs of the saddle-point equations by mathe- 
MATICA. We set d = 5, and take L = 17, T = 50, A = 12n^ , 



and (? = IStt . The ordinate is p 
and the abscissa is . 



AbiM^)/[l-2gb{A'P)], 



where M'^* > A/^*. Here, we compute the free energy 
densities at high temperature. Substituting the saddle 
point iXax = ni^ = M'^ — Q and Qa/Sx = Q into the 
replica partition function H7U|) and the action H71|) . we 
have the following expression for the free energy density 
/ {= F/V) in the large N limit: 



■' 2(3 



In 



- Jd 



25 



a{AP) + j-^a{M-') 



r2\2 



(87) 



As the temperature becomes higher, the intersections 
and become as follows: 



' 



Mi 



T, 



(89) 



where Mq is given by solving the equation 1 = 2gh{M^). 
Thus, the free energy densities at high temperatures in 
both solutions are 



I- 



1 

2p 



M2) 



fiMl* ~ T) 

_^[l + ln(27r)]- — 
2 ^ ^ 2kT 



2/3 



A- 



a{Ml 
9 



2kT 



Jd 
(90) 

(91) 



The free energy density /_ is lower than /+. Perform- 
ing the high temperature expansion without the replica 
method, however, we find that the result is consistent 
with /+ in the leading order. Details of the calculation of 
the free energy density at the high temperature without 
the replica method are relegated to Appendix 1X1 Thus, 



^(pure) 




FIG. 3: The phase diagram. The equation for the boundary 
surface is {l-T^/T^^'^'"'''){l-2gbo) = Ac/Ar"°^ The region 
containing the origin is the ferromagnetic phase. The other 
is paramagnetic. 



the solution {M'^,p*_) should be excluded. This choice 
of the solution (Af^* , is consistent also with the result 
obtained by the functional renormalization group analy- 
sis in the large N limit at zero temperature, as discussed 
in the final section. We also should note that the saddle 
point (A<f^*,p!j_) exists in the region 



1 - 2gb{M^) > 0. 



(92) 



Near the critical point, M becomes small, and then 
the field theoretical description is considered to be appli- 
cable. The integrals (|80|l and (|81|l can be expanded in 
4 < d < 6 as follows: 



a{M^) 
b{M^) 



ao 



- aiA/T^ 
biM''- 



(93) 
(94) 



where oq, ai, 6o, and bi are the same positive constants 
as those of Eqs. (|20|) and 121|) . Inserting the above ex- 
pansions into the saddle point equations and 
we get 

1 = kTaa + (A + 2g)bo - 2g{kT)ao(bo - biM"^-^). (95) 

Putting A/ = 0, we can get the critical line between fer- 
romagnetic and paramagnetic phases: 



fcTcCo + (Ac - 
(1 - 2560) (l 



25)60 - 2g{kTc)aobo = 1, (96) 



(pu 



, (T=0) ■ 



The phase diagram is depicted in Fig. 13 We find that the 
ferromagnetic region is smaller than that in the absence 
of the random anisotropy term. As the strength of the 
random anisotropy increases, the ferromagnetic region 
becomes small. M is rewritten by using Tc and A^ as 
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follows: 
M 



A + 2.9 1 



T 



(pure) 



^l/(d-^4) 



Putting A = Ac, we have 

(l-2g6o)(T-Tc) 



l/(d-4) 



(pure) 



l/(rf-4) 



.(98) 



(99) 



Putting T ~ Tc, we have 

(l-2,g6o)(A- A,) 



M = 



[Ac + (l-2.g6o)(A-A,)]Ai^-°)6i 



l/(d-4) 



(100) 



Thus, the exponent v of the correlation length is = 1/2. 



B. Stability of replica-symmetric saddle point 



We put 



iXax = (M^ - Q) + ifo 

Qafix — Q + Vaflx , 
Xa/3a; = Xa/3 + (^Xa/Sa;- 



(101) 
(102) 
(103) 



In the same way as in the previous section, we expand 
the effective action S'eff up to the second order of SxajSx- 
To study the stability of the saddle point against the off- 
diagonal fluctuations 77^^^ , we calculate the eigenvalue of 
the following Hessian: 



G. 



(q/3)(75) 



Sr]af3kSri^s~k 
Putting the following ansatz (replicon subspace): 



(104) 



(105) 



7=1 



we get the eigenvalue as 



1 



Afe - - 1 - 25 



5V J,,[_.,.].(2.)^^°"-^°" )^^°') 



in the thermodynamic limit. Details of the calculation 
are shown in Appendix ^ If the eigenvalue Afe is positive 
for all k, the saddle point is then stable against the off- 
diagonal fluctuations rjajsk- First of all, putting k ~ 0, 
we can easily investigate the eigenvalue 

^ _ i/^^_2 f 1 

= -^[l-2gb{M% (107) 



The condition that the eigenvalue is positive is given by 



1 - 2gb{M^) > 0. 



(108) 



This is in agreement with the region (|92|l where the sad- 
dle point exists. Thus, the eigenvalue is positive in the 
region of the critical point and over. This result indicates 
that the replica-symmetric saddle point is stable against 
the fluctuation that is induced by introducing the second- 
rank random anisotropy, and therefore it is possible to 
integrate out the fluctuations f]a0x ■ Even though we cal- 
culate the higher order corrections inl/N expansion, we 
cannot find the instability of the replica-symmetric sad- 
dle point against the fluctuation. Therefore, the dimen- 
sional reduction holds for sufficiently large N. 



IV. FUNCTIONAL RENORMALIZATION 
GROUP FOR LARGE N MODELS 



We compare our results with the functional 
renormalization group (FRG) study at the zero 
temperature>i2ii^*i& We search for a consistent FRG 
solution with the 1/N expansion. Details of the analysis 
are given in Appendix [CI In general, a replicated 
Hamiltonian can be written as 

X a 

-yEE^(^--"-^-./5), (109) 

where the function R{z) represents general anisotropy. 
Our Hamiltonian H62|l corresponds to choosing 



R{z) = Az + gz^ 



(110) 



First, we discuss the solutions in the large TV limit. If 
one takes the large N limit, one finds exact solutions of 
all fixed points. We can analyze their stability by solv- 
ing the eigenvalue equation of the infinitesimal deviation 
from the fixed-point solutions. This method is discussed 
by Balents and Fisher^^ for random media. The one- 
loop beta function for a general R{z) has both analytic 
and nonanalytic fixed points. Following the method of 
Le Doussal and Wiese^i^ we find one-parameter family 
of nonanalytic fixed points with a cusp. We obtain an 
asymptotic form of the solution near z = 1, 



R'{z) - i?'(l) + ^2R'(l)[e/A- R'{1)]{1- z). 

Our analysis shows that all physical nonanalytic fixed 
points satisfying the Schwartz-Soffer inequality^^ have 
many unstable modes. 

In addition to the nonanalytic fixed points, we find 
four analytic fixed points given in Eq. H110|l with 
(A,g) = (0,0),(e/A0),(0,e/(2A)), and (-e/A, e/(2A)), 
where e = d — 4. The last one is unphysical since A < 0. 
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The triangle defined by the other three fixed points corre- 
sponds to ferromagnetic region on the T ~ plane in Fig. 
El In fact, the vertex (A,g) = (aF"°\ 0) = (l/6o,0) 
corresponds to H97|) with Tc = 0. It is easily seen that 



1 

b~o 



A 



(111) 



where A is a momentum cutoff. Thus the dimension- 
less quantity (A,g)A'^ is equal to the analytic fixed point 
{e/A, 0). Similarly, the vertex (0, l/(26o)) corresponds to 
the fixed point (0,e/(2^)). Therefore, we find that the 
phase diagram at T = obtained by the large N limit 
is understood by the functional renormalization group 
method. Furthermore, the stability analysis in Appendix 
Elshows that {e/A, 0) is singly unstable, where the unsta- 
ble mode corresponds to deformation along A axis. The 
origin (0, 0) is fully stable while (0, e/ {2A)) is fully unsta- 
ble. Therefore, in the large N limit, a phase transition 
at T = is governed by the singly unstable fixed point 
(e/j4, 0) yielding dimensional reduction. The correspond- 
ing flow in the two dimensional coupling constant space 
is depicted in Fig. 0] 




FIG. 4: The renormalization group flow for the couplings A 
and g in the large N limit. 



Next, we discuss the model with a finite N. We con- 
clude that the dimensional reduction holds for sufficiently 
large N. We study the singly unstable analytic fixed 
point found in the large N limit. By the discussion in 
Appendix [Q the fixed point to control the phase transi- 
tion has 



i?'(l) = 



A{N - 2) 



or 



0. 



At this stage, we find only two possibilities. The expo- 
nents of the correlation function becomes those given by 
the dimensional reduction 



rj — 1] 



N-2 



or the trivial ones 



77 = 77 = 0. 



(112) 



(113) 



We obtain the subleading correction to this fixed point 
solution 



m = -J 



1 1/2 

2 + 2^(^ 



2z] 



(114) 



We analyze the stability of this analytic fixed point yield- 
ing the dimensional reduction by solving the eigenvalue 
equation of the linearized beta function in Appendix 
mi There are many unphysical modes which diverge 
in the interval — 1 < z < 1. These are not generated 
in the fiow, thus we eliminate these unphysical modes 
by choosing the integral constants. These unphysical 
modes correspond to the infinitely many relevant modes 
pointed out by Fisher. ^■^ In our solution of the eigen- 
value equation, we find the same eigenvalues calculated 
by Fisher up to the order 1/A'^, if we correct an expression 
given there by adding an overlooked term. We discuss 
this problem in Appendix |0 The analytic fixed point 
(e/^,0) has slightly relevant operators with dimension 
less than 2/A^, which give deformation of the coupling 
5R'{z) ~ (1 - z)-" with < -a < 2/N. Here, we dis- 
cuss this subtle problem of the slightly relevant operators. 
First, we assume that the initial coupling constant i?"(l) 
in the renormalization group equation is finite. In this 
case, this fixed point behaves as a singly unstable fixed 
point in the following reason. By Fisher's representation 
of the renormalization group equation, -R'(l) and -R"(l) 
satisfy 

dtR!{l)^{^ - d)R'{l) + A{N - 2)i?'(l)2, 
dtR"{l) = {^ - d)R"{l) + A(6i?'(l)ii"(l) 

+ {N + 7)R"{lf + R'{lf). (115) 

For a small initial value of i?"(l) the flow of i?"(l) stays 
in a compact area. The flow in the two-dimensional cou- 
pling constant space is qualitatively the same as in Fig. 

H 

In this case, if i?'(l) takes a critical value by tuning 
the coupling constant A or the coupling R{z) flows to- 
ward the analytic fixed point with a finite i?"(l). Then, 
the fiow does not generate the relevant mode with an ex- 
ponent < — a < 2/N from an initial function with a 
finite -R"(l). This analytic fixed point controls the phase 
transition, and therefore the critical behavior obeys the 
dimensional reduction. Since this analytic fixed point 
exists for > 18 as pointed out by Fisherji^ the di- 
mensional reduction occurs for A^ > 18. In this case, 
the critical exponents of correlation function are given 
by (|112|l . This result agrees with our simple 1/A^ expan- 
sion. Next, we consider that the initial coupling constant 
R"{1) is not finite. We assume 

R'{z) = C(l -z)-", 

with < —a < 1. Since R"{1) diverges, already at 
the initial stage the coupling constants are infinitely far 
from the analytic fixed points for any small C. We can- 
not justify whether or not the continuum field theory 
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approximation induces such a mode in the initial func- 
tion. TheoreticaUy, however, we can consider such a 
model. The renormalization group transformation can 
generate a term proportional to (1 — Since 
—2a — 1 < —a, the successive transformation may pro- 
duce less power. Eventually, the flow generates a relevant 
mode with — a < 2/N and the flow cannot reach the ana- 
lytic fixed point by tuning the parameter A and g. Since 
all fixed points are unstable except the trivial one, the 
flow reaches the trivial fixed point directly in the mass- 
less phase. In this case, we obtain only trivial critical 
exponents H113|) . This second possibility does not agree 
with our 1/A^-expansion method. Therefore, only con- 
sistent result with the 1 /N expansion is the dimensional 
reduction. 



V. DISCUSSION AND SUMMARY 

In this paper, we have studied the random field 0{N) 
spin model including the second-rank random anisotropy 
term. We have studied the effect of the second-rank ran- 
dom anisotropy on the critical phenomena of the ran- 
dom field O(A^) spin model in 4 < d < 6, by use of the 
replica method and the 1/iV-expansion method. The off- 
diagonal fluctuations are induced through the Hubbard- 
Stratonovich transformation for the second-rank random 
anisotropy. We have computed the saddle point under 
the assumption of the replica symmetry, and have stud- 
ied the stability of the replica symmetric saddle point 
against the off-diagonal fluctuations which are induced 
by the second-rank random anisotropy. Our criterion to 
judge the stability of the system is identical to the stan- 
dard one used by de Almeida and Thoulessii^ It is based 
on the stability of the saddle point of the auxiliary field 
introduced to calculate the partition function explicitly. 
We find that the eigenvalues of the Hessian around the 
replica symmetric saddle point are positive definite, and 
thus the Gaussian integration over auxiliary field can be 
performed. The instability is not observed in higher or- 
der correction in the 1/N expansion. Consequently, we 
conclude that the replica-symmetric saddle point is sta- 
ble for a second-rank random anisotropy with the order 
1/N and the dimensional reduction holds for sufficiently 
large N. 

This result is inconsistent with that obtained by 
Mezard and YoungiiSi Since the SCSA equation gives 
the precise two-point correlation function up to order 
1/7V, their replica-symmetric two-point correlation func- 
tion agrees with ours. Nonetheless, they conclude that 
the replica-symmetric correlation function is unstable 
against a deviation of the correlation function by treating 
the free energy as a functional of two point correlation 
function. Their criterion for stability differs from the de 
Almeida-Thouless one, although it looks the same. They 
optimize the free energy by choosing the two-point cor- 
relation function freely. On the other hand, in our anal- 
ysis, a two-point correlation function can be deformed 



only through changing a saddle point of the auxiliary 
field, and then it cannot be deformed freely. This is the 
essential difference between two theories. We consider 
either that the instability shown by Mezard and Young — 
is just apparent, or that their method includes some non- 
perturbative effects other than the 1/iV expansion. For 
the latter possibility, we should justify that the free en- 
ergy can be optimized by a correlation function with no 
constraint. 

We have checked the consistency between the large N 
analysis and the renormalization group fiow by showing 
that the phase boundaries obtained in those methods are 
consistent in 4-1- e dimensions. As pointed by Feldmanji^ 
the critical phenomena near the lower critical dimension 
is governed by the nonanalytic fixed point by the ap- 
pearance of the cusp, and then the dimensional reduc- 
tion breaks down for some small N. For large N, how- 
ever, we show that the functional renormalization group 
method studied by Feldman allows us to perform the 
1/iV expansion. We find all fixed points which consist 
of analytic and nonanalytic ones in the large N limit. 
On the other hand for N < 18, it is known that there 
are no nontrivial analytic fixed pointsii^ By solving the 
eigenvalue problem for the infinitesimal deviation from 
the fixed point, we find that the nonanalytic fixed points 
are fully unstable. We search for consistent solutions of 
the renormalization group with the 1/N expansion. If 
the initial R"{1) is finite, the nonanalytic relevant modes 
cannot be generated. In this case, the unique analytic 
fixed point practically behaves as a singly unstable fixed 
point, which gives the dimensional reduction. This result 
agrees with the stability of the replica-symmetric saddle- 
point solution in the 1/iV expansion. Thus, we conclude 
dimensional reduction occurs. 

Our result also agrees with a recent study of the ran- 
dom field 0{N) model by Tarjus and Tissier. They study 
the model by a nonperturbative functional renormaliza- 
tion groupjSSi Although their work to obtain a full so- 
lution is in progress, they give a global picture in a d-N 
phase diagram and discuss the consistency of their results 
with those by some perturbative results. They propose 
a scheme to fix a phase boundary of the phase where the 
dimensional reduction breaks down. Using an approxi- 
mation method, they show that the phase is in a compact 
area on the d-N plane. 

Here, we comment on the model in dimension less than 
4. The l/7V-expansion method shows that the model has 
a massive paramagnetic phase only. Also, the functional 
renormalization group method for negative e — d — A 
shows that there are no nontrivial analytic fixed points. 
The trivial fixed point and nonanalytic fixed points are 
unstable for c? < 4. Our large N analysis indicates that 
the nonanalytic fixed points are unstable, and therefore 
only a massive phase exists. This result agrees with Feld- 
man's resulti^ that the correlation length is finite always 
for N > 10. 

Finally, we comment on the critical behavior near the 
upper critical dimensions. In a recent work,^^ the dimen- 
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sional reduction has been shown by a perturbative renor- 
mahzation group in a couphng constant space near the 
upper critical dimension in the random field Ising model 
at one-loop order. This result is also consistent with that 
obtained by Tarjus and Tissier,^'' This study can be ex- 
tended to the 0{N) model and the result agrees with the 
1/A^ expansion. These studies suggest that the large N 
limit may be applicable to the model with a small N near 
the upper critical dimensions. However, it is a nontrivial 
problem whether or not the dimensional reduction holds 
near the upper critical dimension for a small N. Further 
studies are needed. 



Performing the calculation of the measure J DS, we have 

y = n 2r(jv + i) - exp[y.(iv)], (A6) 



s{N) ^[l + ln(27r)], 



(A7) 



for > 1. 

We study the behavior of the free energy at high tem- 
peratures. We expand the partition function in /3H up 
to the second order: 
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APPENDIX A: FREE ENERGY AT HIGH 
TEMPERATURE WITHOUT REPLICA METHOD 



Z = / VSe-^" 



e^^^^^l-im + ^imf)], (A8) 



where the angular brackets (• • • ) stand for 



_ JVSfiS) 



J VSfiS). (A9) 



The Hamiltonian is given by 



Then, In Z is 



H = —J ^ Sx ■ Sy — ^ hx ■ Sx — ^(^K • Sx)'^ 

{x,y) X X 

N N 

(x.y) i=l 



X i—1 



N 



/ , / , "£C "-a; '^x "-"a; 



X l.J 



N 



N 



{x,y) i=l X i=l 



N 



/ , / , "'2.x '-'x '-'x ■ 



(Al) 



Here, h!i^ and h2^^ are the random field and the second- 
rank random anisotropy, respectively: 

= 0, \h^h^'^]^M,,5xy, (A2) 
= 0, = \{5,u5,i+5u5,u)5xy. (A3) 

The partition function is 



^ = (n / dSxSiSx'' - 1) 



,-0H 



Here we put 



n r dSxSisi ~i)= fvs 



(A4) 



(A5) 



InZ ~ Vs{N)^ {(3H) + -{I3H;I3H), (AlO) 



where 



{i3H-i3H) = {mf)-mf 



(All) 



Using the identity (Sx ^Sy^) = SxySij, and Eqs. HA2p and 
^311. we have 



1 



f = -fr.Z] 



NkT,^ , M J d NA + Ng 

~ 1 + ln(27r -. 

2 ^ ^ 2NkT 2NkT 

(A12) 

According to the redefinition of the parameters (|59|l . the 
free energy density is rewritten as 

/ = -ll^^z] 

P 

kT, , , J^d A + g , 

~ l-hln(27r —. (A13) 

2 ^ ^ 2kT 2kT ^ ' 

This is in agreement with Eq. (|91|) . 



APPENDIX B: DERI VATI ON OF EIGENVALUE 

(IToet 

In this appendix, we give the details of the calculation 
of the eigenvalue (|105|l . 
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The second-order term of Sxafix and rjapx for the ef- 
fective action S^s becomes 



Tr- 



x< X 



1 



-Sx — - — 



-5x 



-f f d'xY^vl^x- 



(Bl) 



In the momentum representation, the second-order cor- 
rection of the action 5cff is rewritten as 



N 



c« — -r 



fce[-x,7r]"i (271-)'' Vf- 



+ 2 ^ ^ (-ak'niS~kGa(^S) 
a—1 7<(5 



Val3kV'f5,-kG(^al3){jS) j ■ (B2) 



a<f3j<5 



Here 



(B3) 
(B4) 



G, 



-t-Lrg * (JQ 

EC (a ^5), 
i 

'2 



G^ * 



(Go'' * G^" -f G^" * Go^ 



'0 



G 



(q/3)(q/3)^ 



GT * Gf) 



1 



(G^ * G^ 



+G(^'^ * G^' 



^0 



^0 



Gf) 



S0 



(a/3)(a5)--(Go * Gg + Gq * Gq 

+ Gf * Gf^" + Gf * Go^) 

G(a;3)(7<5) = -(Go'' * Go" -I- Gq'^ * Go'' 



(B5) 



(B6) 



(B7) 



(B8) 



-Gf * Gl"" + G:?-^ - ^''"^ 



E7^ (a^/3,7,,J,/?^7,57^5). (B9) 



Ga/3 y^7(5 
* ^-'o 



(BIO) 



We shall find the eigenvalues of the matrix G (Hessian) 



Q — i {^07} {Gq,(^5)} 
{G(a/3)7} {G(q,^)(^5)} 



(Bll) 



Let /2 be the eigenvector associated with the eigenvalue 
A; 



where t]-^^ — 'HS'y- Applying G to ^, we obtain 
(G/i)a = Aeo,+Bj2'^-r+<^Yl 

7:/a 77^<^ 



(B12) 



V 



(B13) 



)+V £j + 'Pv 



-2Q Y^^^^ + ^Y. E 



n n 



7^a,/3 



7^Q,/3 (5^7, a, /3 



(B14) 



To find the solution of Eqs. (|B13() and HB14|I . we use the 
following ansatz according to Ref. 10: 



6q — ^ ^ ^07 — 

7=1 



(B15) 



for all a. Under this ansatz, Eq. (|B14(I gives a nontrivial 
solution 

{Gtl)a.t3 = {r-2Q + n)i^^p - Xii^p. (B16) 
Therefore, we get the eigenvalue (|106|) . 



APPENDIX C: FUNCTIONAL 
RENORMALIZATION GROUP STUDY FOR 
CRITICAL PHENOMENA OF RANDOM FIELD 
0(iV) SPIN MODEL IN 4 + e DIMENSIONS 

In this appendix we study the one-loop beta function 
derived by Fisheni^ for a general random disorder R{z) 
at zero temperature: 

dtR{z) = (4 - d)R{z) + A ^2(7V - 2)R'{l)R{z) 

-{N - l)zR'{l)R'{z) + (1 - z^)R'{l)R"{z) 
+ - 2 + _ R'{z)R"{z)z{l - z^) 



+ ^[R"iz)ril-zY 



(CI) 



Here, i = In^ with I being the length scale specifying the 
FRG and A = Sd/{2Try,P. 



14 



1. General properties of fixed points 



2. Large- limit 



The fixed point condition of the renormalization group 
determines properties of the function R{z). Here we dis- 
cuss possible asymptotic behaviors of R{z) near z = 1. 
The first derivative of the fixed point equation with re- 
spect to z is 

[(4 ~d)/A+{N- l)]i?'(l)i?'(z) + zR'{zf 

-{N + l)zR'{l)R"{z) + iN-2 + Sz^)R'{z)R"{z) 

+ (1 - z^)R'{l)R"'{z) - z{l - z^)R'{z)R"'iz) 

-3z(l - z^)R"{zf + (1 - zyR!'{z)R"'{z) 

-{I- z^)R'{z)R"{z) = {). (C2) 

If we assume asymptotic behavior of R'{z) near z — \, 

R\z) = R'{l) + C{l-z)-< + ■■■ , (C3) 

with < 7. To discuss a cuspy behavior of R{z) at 
z = 1, we consider only 7 < 1. The condition ljC2|) gives 
the following constraint: 

[(4-d)M+(Ar-2)i?'(l)]i?'(l) 

-(7^7(472 + 47 + A/- - 1)(1 - zf"'-^ = 0. (C4) 

For 7 7^ 1/2, this constraint gives 
1 



7= -(-1 + V2-iV) or C = 0, 



and also 



R'{\) 



d-A 
A{N - 2) 



or R'{\) = 0. 



Here, the former case shows the dimensional reduc- 
tion. The formulas for the critical exponents obtained 
by Feldmanji^ 

r; = Ai?'(l), f] ^ A{N - l)R! {1) - e, (C5) 

enable us to obtain 

d-A 



N -2 



(C6) 



In this case, no 7 is allowed for any N. For 7 = 1/2, 
the parameter R'{1) can change continuously depending 
on the constant C. Therefore, only 7 = 1/2 allows di- 
vergent R"{1). Only this case does the nontrivial critical 
behavior differ from the dimensional reduction. Since the 
initial value R{z) of the renormalization group equation 
(|C1|) is an analytic function, the flow of R"{1) should 
diverge for the breakdown of the dimensional reduction. 

The same discussion for z — ~1 can be done. The only 
possible singularity is 

R'{z) = R'{-l) + C{l + zy/^ + --- . 

If C = 0, then we have 

R'i-l) = iN-'S)R\l)-^-^ or i?'(-l) = 0. 



In order to take the large N limit, we multiply both 
sides by N and rescale NR R. The beta function 
becomes 

dtR{z) = {A - d)R{z) + a{2R'{1)R{z) - zR'{l)R'{z) 
1 



+ -[R'{z)r]+Oil/N). 



3. Fixed points 



(C7) 



Following the method given by Balents and Fisher^ii 
we consider the flow equation for R'{z) instead of that 
for R{z). Taking the derivative with respect to z and 
introducing u{z) defined by 



R'{z) = -u(z 
A ^ 



(C8) 



the fixed point equation for (|C7p becomes 

(a - l)u{z) - zau'{z) + u{z)u'{z) = (C9) 

in the large N limit. Here we define a — u{l). First we 
solve it when a — 1. In this case, u{z) satisfies u'{z) = 
or u{z) = z. If u'{z) = then u{z) = 1 since a — u{l) — 
1. Thus 



(CIO) 



where the constant term —e/{2A) is determined by ljC7p . 
On the other hand, in the case of u{z) = z, 



(Cll) 



Next, we turn to the case of a = 0, where u{z) satisfies 
u{z) = or u'{z) = 1. The former case is R{z) = 0, which 
corresponds to the pure theory. The latter becomes 



namely, 



u{z) = (z - 1); 



(C12) 



(C13) 



Those analytic fixed points were first obtained by 
Feldman.^^ 

Next we consider a general case. If a 7^ 0, 1, 



du {a — l)u 
dz za — u 



(C14) 



Taking the inversion we regard z as a function of uJ^ 
One gets 



dz a z 1 



du a— 1m a— 1' 



(C15) 
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FIG. 5: A schematic graph of u{z). Since the derivative of u 
is ill defined on u = az, the solution terminates on this line. 
The above graph represents two solutions meeting at (1, a). 



which is easily integrated The result is 



(C16) 



where C is a constant. Since z{u) satisfies z{a) — 1, C is 
determined uniquely as 



z(u) = (1 — a) 



a/(a-l) 



(C17) 



Now we revert (|C17p to the solution u{z) for (jC14p . Be- 
cause z{u) takes the maximum value 1 at m = a, u{z) 
is double valued as we show in Fig. |S| It is seen from 
(|C14|I that du/dz is ill defined on m = az. Therefore the 
lower branch terminates at the origin, so that it should 
be continued to the region — 1 < z < 0. This is possible 
only if a/{a — 1) is a positive integer. 

Nonanalytic behavior near z = 1 is clarified as follows. 
Set 



u = a + Su 
and assume that \Su\ < \a\. From !KJT7^ . 



i{z) = a ± A/2a(a - 1)(1 - z) H 



(C18) 



(C19) 



Note that the plus (minus) sign in front of the square 
root means to take the upper (lower) branch. Nonana- 
lytic behavior is seen at z = 1. Since the function u(z) 
should be real in |z| < 1, a satisfies a (a — 1) > 0. Further- 
more, fj — ea/N should be nonnegative due to physical 
requirements [see (|C5|l ]: hence 



a > 1. 



4. Stability of the fixed points 



(C20) 



Next we investigate the stability of the solutions. Let 
u* be a fixed point solution: 



l)u* 



zau + u u 



0. 



(C21) 



To study the stability of u* , let u ^ u* + v. Inserting 
this into IIC9|l and keeping up to the linear terms of i;, 
we consider the following eigenvalue problem: 

[{a - 1) + u']v + {u~ za)v' + (u - z-it>(l) = \v. (C22) 



Here we omit the asterisk from u* for brevity. Normal- 
izing V appropriately, we can take viV) = or v(\) — 1. 
We begin with the analytic cases. 



R{z) = e{z - 1/2) /A 



In this case, a — 1 and u{z) = 1. Then (jC22|l becomes 



(1 - z)v' + b = Xv, 



(C23) 



where b represents v(l) taking or 1. When b = the 
solution is 



v{z)^C{l-z) 



-A 



(C24) 



where A < because of the initial condition b = v{l) = 0. 
On the other hand, when & = 1, a general solution is 



/ A / A-i+c(l-z)-^ 
^(")= ln|l-z 



(C25) 



(A ^ 0), 
(A = 0). 

Here the condition 6=1 requires that A = 1 and c = 
0. In conclusion, the allowed value of A is A < or 
A = 1. This shows that the fixed point solution is singly 
unstable. 



b. R{z) = ez^/{2A) 

In this case, a = 1 and u = z; hence (|C22|I is simplified 
to V = Xv for 5 = 0, 1. It means that A = 1 for every 
deformation, so that this fixed point is fully unstable. 
This is also true for any finite TV. 

c. R = 

Since a = and u = in this case, ljC22l) is —v — Xv, 
which means A = — 1 for any v; thus the trivial fixed 
point is fully stable. 

d. R{z) = e{z - 1)V(24) 
Here, a — and u — z — 1. The eigenvalue equation is 



(C26) 



z)v' ~b = Xv, 

which can be solved in a similar way as for (|C23|I . The 
result is 



for = and 



v{z) 



v{z) = C{l-zy 



-A-i + c(l-z)^ (A 7^0), 
ln|l-z| (A = 0) 



(C27) 
(C28) 



for 6 = 1. Therefore the allowed values of A are 

A > and A 1. (C29) 
Therefore it is unstable. 
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e. Nonanalytic case 

Next we proceed to the nonanalytic case. Using (|C14p . 
we regard v as function of u. Then (|C22p is written as 



dv 

— + f{u)v + g{u)b^Q, 
du 



(C30) 



where 



/(«) 
9{u) 



A ^\ 1 1 



a— 1 J u az — u 

z — u 



(1 — a){za — u) 
A general solution of ljC30|) is 



where 



(6 = 0), 



F{u) = / f{u)du 



(C31) 



(C32) 



(C33) 



Let us compute F{u). Since z is given as a function of u 
by ljC17p . we can write 



du 



du ( 1 



1 - a \u^l'y°'-^') - 1 



where 



,(C34) 



(C35) 



Thus, using the ambiguity of the constant term of F{u), 
we get 



A — a 



F{u) = Inu + In |1 - u^^^''~^^\. (C36) 

a — 1 



Therefore, 



e-Hu) ^ ^(a-A)/(a-l)|^ _ ^il/(a-l)|-l. (C37) 



When 6 = 0, w is proportional to (|C37|) . which becomes 
singular at u = a, i.e., z = 1. Hence, there are no non- 
trivial solutions satisfying 6 = 0. 

Next we consider the case 6=1. From (|C31|I and 
(ESU), we get 



e^^-^giu) = ±- 



~A/(a-l)-l 

(1 — a)a 



(C38) 



Note that the plus sign is taken for the upper branch and 
the minus for the lower branch. Inserting this into HC32|I . 
we get 



v{u) = ' ^ 



(l-a)(l-ui/(a-i)) — "^ 



Here the constant terms are chosen to satisfy v{u{z)) 1 
as z —^ 1. Thus, the deviation v{u) from the upper 
branch is finite for any A, because it > 1. On the con- 
trary, v{u) from the lower branch may diverge at u = 
and —1. We need a constraint on A for v{u) to be fi- 
nite. We find that the lower branch with a = 3/2 can 
be extended to — 1 < < 0, and that v{u) remains fi- 
nite for A = 1 or negative integers; namely, the lower 
branch with a = 3/2 is singly unstable. However, this 
fixed point solution is unphysical because it does not 
satisfy the Schwartz-Soffcr inequality 2ri > fjm^ This 
inequality requires a = 1 + 0{1/N). Other physical 
lower-branch fixed points satisfying the Schwartz-Soffer 
inequality have many relevant modes of 0{N). 



5. Subleading corrections 

a. The stable fixed point and critical exponents 

Here, we calculate the subleading correction to the an- 
alytic fixed point R{z) = (e/A){z — 1/2) and the eigen- 
functions. We expand the fixed point solution 

R{z)^^Ro{z) + ^R,{z) + o(^^y (C39) 

and calculate the subleading correction Ri(z). Substi- 
tuting this expansion into ljCl|) . we obtain 

dtRi (z) = (4 - d)Ri (z) + A (^2R[{l)Ro{z) (C40) 

+2R'oil)Ri{z) ~ zR[{l)R'o{z) 
-zR'oil)R[{z) + R[{z)R'o{z) 
-AR'oil)Roiz) + zR'o{l)R'o{z) 

+ il-z^)R',{l)R'^{z) + liz^-2)R',izr 

-R',{z)R'^{z)z{l - z') 

+ i<(z)2(l-z2)2 

We substitute the unique singly unstable fixed point so- 
lution 



Ro{z)^j(z^l 



into the above equation; then we obtain a fixed point 
equation for the corresponding correction Ri{z), 



(1 - z)R[{z) + Ri{z) + (1 - z)R[{l) 

(C41) 

We obtain the following unique solution of this equation: 



+ l(lz^-3z+l)=0. 



Ri{z)^—{z' + 2z). 



(C42) 

Fisher indicated that this fixed point exists for N > 18. 
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b. Stability of the analytic fixed point 

We substitute the analytic fixed point expanded in 1/iV 
into the eigenvalue equation for an infinitesimal deforma- 
tion of the coupling function 

(1 - zf{l + z)v"{z) + {l-z){N -Az- 2)v'{z) 
+ {2z^ N\)v{z) + {N -2)v{l) = Q. (C43) 

First, we study the equation for v{l) = 0. Solutions of 
this equation have regular singular points z = 1 and —1 
for the interval — 1 < ^ < 1. Therefore, we can obtain the 
solutions in the following expansion forms around z = 1 : 



i;(z) = (l-z)-"^a„(l-z)", (C44) 

n=0 

and around z = — 1 

oo 

= (l + z)'5;^6„(l + z)". (C45) 

Substituting these forms into the eigenvalue equation, we 
require that the coefficient of the lowest order vanishes. 
This requirement gives the indicial equations for the ex- 
ponents a and /3 

2a^ (iV - 4)a + 2 - iVA = 0, /3(2/3 + iV) = 0, (C46) 

which have solutions 



A- N ±^/W^m + W\ ^ N ^ 
a± = ^ , (3^- — ,0. 

(C47) 

The coefficient of an arbitrary order satisfies the following 
recursion relation: 

2k{k — a± + Q!ip)a^ — (a± — k){a± — k — \)a^_^ = 0, 

for fc = 1, 2, 3, . . . . By solving this recursion relation, the 
expanded solution can be written in the Gaussian hyper- 
geometric function as follows: 



/ iV 1 - z 

= FI l-a±,2-a±,3-2a± - y;^— 

(C48) 

Solutions with a > or /? < diverge at z = 1 or —1, 
and they are unphysical. To obtain a finite solution for 
the interval — 1 < z < 1, we construct a general solution 
as a linear combination of two solutions, 

oo 

^;(z)=C+(l-z)-"+Ea+(l-^r 

oo 

+C_(l-z)-"-Ea-(l-z)". (C49) 



We can eliminate the divergent solution with j3 = ~N/2 
at z = —1 by choosing C± for a requirement 1)| < 
oo. Also the finiteness of v{l) requires a± < 0, then we 
obtain a condition on the eigenvalue 



(C50) 



This condition on A implies the existence of slightly rel- 
evant modes at this analytic fixed point. In addition to 
these modes, we find one relevant mode for ^(l) 
with A = 1 by solving the eigenvalue equation, as well 
as in the large N limit. This fixed point yielding dimen- 
sional reduction seems to be unstable except in the large 
N limit. There is no singly unstable fixed point gener- 
ally. The only stable fixed point is the trivial fixed point. 
In a limited coupling constant space where R"{1) is fi- 
nite, however, the analytic fixed point is singly unstable. 
Then, dimensional reduction occurs in such models with 
a finite R"{1) as initial coupling constant, as discussed 
in Sec. |IV| 

Here we comment on the infinitely many relevant 
modes pointed out by Fisher They are included in the 
following series in our solution (|(y49p : 



a_ l-/c, (A: = 3,4,5,...) and C+ = 0. 
These belong to the eigenvalues 

which are positive for sufficiently large k. These agree 
with the eigenvalues obtained by Fisher, although we 
should add a term 2nkP2Pk missed in Eq. (C6) of his 
paper. Since these relevant modes diverge at z = — 1, we 
have eliminated them as unphysical modes, as discussed 
above. 

APPENDIX D: INTEGRALS 

We restrict ourselves to 4 < d < 6. 

, . d'^k 1 



d'^k 1 

d'^k 1 



.(Dl) 



We put 



d'^k 1 



el-TT,^]-^ {'^'^r -JAfc 



ao > 0. (D2) 



We calculate the second term. Putting k mk/\/j, 
and using the approximation — A^^j,,^ ~ m'^k'^/J for 
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<C 1, we have 



m k I J for m <C 1, we have 



kel-^^M-" (2'^)'' (-JAfe)(-JAfc +m2) 



{2-K)dJd/2 

(27r)''j2 d - 4 



■K\fJlm k-d-S 

dk- 



P + l 



(D3) 



where Sd ^ 27r'^/2/r(d/2). Thus, we have the following 
expression for a{m'^): 



d'^k 



(D4) 



d'^fc 1 



i^Trf JAfe)(- JAfe + m2)2 



fee 



[-.,.]^ (2^)'^ (-JAfc)2(-JAfc +m2)2 



We put 



d'^k 1 



(D5) 



= 6o > 0. (D6) 



We calculate the second and the third terms. Putting 
k — !■ nik/^fj , and using the approximation —A^^^,^^ ~ 



2m^ 



1 



•^-4 ' dk- 



OQ, fTT\/j/m ud-3 

—m'^-'^ / dk 







(2,^)^7^/2 



and 



m 



d'^k 



ke[-^M' ^"^^y {-JAk)^-JAk + m2)2 

Jd d-i I 



{2n)dJd/2 



dk- 



(fc2 + 1)2 



Sd (6" rf)7rcosec(^d/2) ^^^_4 ^ 0{m'^), 



l^2n)djd/2 



(D7) 



Then, 



2m" 



d'^k 



1 



A:G — TT.TT P 

{2Tr)d (_JAfe)(-JAfe 
d''fc 1 
fee[-.,.]^ (27r)'^ (-JAfe)2(-JAfc + m2)2 
5d {d - 2)7rcosec(^d/2) ^ ^^^3^ 



(2Tr)dJd/2 4 



(D8) 



Thus, we have the following expression for b{m^): 

b{m^) ~ bo-bim'^-^. {D9) 
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